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ABSTRACT 


n e q n at vapor pressure 


Two formulas, one for the characteristic time of 
filling a void with the vapor of the surrounding liq- 
uid, and one of filling the void by diffusion of the 
dissolved gas in the liquid, are derived. Based on 
this analysis, it is seen that in an oil film bearing 
operating under dynamic loads, the content of the 
cavitation region should be oil vapor rather than the 
air liberated from solution, if the oil is free of 
entrained air. 

NOMENCLATURE 

C mass concentration of gas 

C s saturation value of C 

D coefficient of binary diffusion 

H Henry's constant 

h thickness of liquid film 

i con d number of molecules condensing on unit surface 
area in uni t time 

Miff mass ® as diffused i nt0 the void from unit 
surface area of liquid in unit time 

Mvap num ^ er °f molecules evaporating from unit sur- 
face area in uni t time 


p pressure 

p s pressure corresponding to C s 

R radius of circular void 

r radial distance in plane polar coordinate sys- 

tem, measured from center of void 

universal gas Constant, «= 8.314 kJ/(kmole • K ) 

T absolute temperature 

t time coordinate 

Miff characteristic time of diffusion, Eq. (22) 

Myn characteristic time of dynamic operation 

Mvap characteristic time of evaporation, Eq. (4) 
v m mean molecular speed 

E dimensionless length coordinate 

p density 

t dimensionless time coordinate 

Subscripts : 


k 

Boltzmann 

constant, - 1 . 380622xl0“ 23 J/K 

a 

thickness 

of diffusion layer 

M 

molecular 

weight 

m 

mass of a 

molecule 

n 

number of 

molecules in unit volume 


g of gas 

i of liquid 

INTRODUCTION 

The cavitation of oil films in bearings operating 
under dynamic loads is not well understood. One of the 
unresolved questions is whether the cavitation bubble 
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contains air, or vapor, or is just an empty void. This 
question is closely tied with the behavior of the bub- 
ble, e.g., whether it collapses or sustains itself dur- 
ing the high pressure portion of the pressure cycles. 
Understanding the bubble behavior is indispensable to 
the construction of a sound theory of dynamically 
loaded bearings. At present, reliable design informa- 
tion is still lacking for important applications of 
dynamically loaded bearings, such as the squeeze film 
damper . 

Air and vapor cavities arise in oil film bearings 
in several ways (Dowson and Taylor, 1979). Air from 
the surrounding atmosphere may be entrained in oil when 
subambient pressure occurs in the bearing clearance; 
the dissolved air in oil may be liberated when the 
pressure falls below the saturation pressure of the 
solution; and oil molecules may evaporate when the 
pressure falls to the vapor pressure of oil. Air 
entrainment is device dependent; it can be dealt with 
only according to specific circumstances. On the other 
hand, liberation of dissolved air and oil evaporation 
are processes dependent only on fluid properties; they 
can be treated in general terms. In the following, the 
latter two processes are analyzed and their respective 
time scales determined. By comparing these time scales 
with that of bearing operation, one can judge the rele- 
vance of these processes with reference to the cavita- 
tion in dynamically loaded bearings. 

CHARACTERISTIC TIME OF EVAPORATION 


_Q_ _ l _ e ' l/t evap 
n eq 

The parameter t e y 3 p signifies the time needed to 
build up the pressure in the void to 63 percent of the 
vapor pressure, and may be taken as the characteristic 
time of evaporation. It is interesting to see from 
Eq. (4) that t eva p is essentially the time a molecule 
takes to travel across the void, and is independent of 
the vapor pressure of the liquid. With the mean molec- 
ular speed given as (Present, 1958, p. 78) 



which is based on the Maxwell distribution, this char- 
acteristic time can be readily calculated. 

In the following, an oil-like substance, squa- 
lene, of known molecular structure (Morales and 
Buckley, 1988) is used for the calculation. The 
needed properties of squalene are given in Table 1 . 

From the observation of the cavitation in a dynami- 
cally loaded bearing (Sun and Brewe, 1990), the cav- 
itation region was large and covered a fraction of the 
bearing surface. If the typical size of the cavita- 
tion bubble is taken to be R = 10 mm and the operating 
temperature to be T * 293 K, then Eqs . ( 6 ) and (4) 
yield v m - 120 m/s and t eva p * 0.167 ms. 


Consider that a void is suddenly created in a liq- 
uid film, the molecules start to evaporate from the 
liquid surface into the void. After some time, the 
pressure in the void is built up to the vapor pressure 
of the liquid, which signifies the equilibrium between 
the flux of evaporation and flux of condensation. The 
vapor pressure of the liquid is assumed low. Hence, 
one may consider (Langmuir, 1913) that the flux of 
evaporation is independent of the presence of vapor in 
the void; and that the flux of condensation is likewise 
not affected by the flux of evaporation, and is deter- 
mined only by the pressure in the void. Since the flux 
of molecules is one-fourth the number density times 
mean molecular speed (Present, 1958, p. 22 ), one may 
write: 


i - - n v , j j - 7 nv ( 1 a , 1 

J evap 4 eq m J cond 4 m v 

Consider that the film has thickness h and the void 
is circular with radius R. Then, the change in the 
number density in the void may be expressed as: 


CHARACTERISTIC TIME OF DIFFUSION 


The solubility of gas in liquid is given by 
Henry’s law: 



. 2 

H 


which states that the mole fraction of the gas in the 
solution is proportional to the pressure of the gas 
outside the solution. The coefficient appearing in 
this relation is called Henry’s constant whose value 
depends on the particular combination of gas and liq- 
uid. Henry's law applies to weak solutions, i.e., 


Hence, Eq. (7) may be written approximately as: 


(8) 


irR 2 h 4? « 2 irRh( j 

dt VJ evap 


J cond J 


Combining Eqs. (l) and (2), one obtains: 


dn n eq n 


dt 


evap 


where 


( 2 ) 


(3) 


C 



Consider a liquid film saturated with a dissolved gas. 
According to Eq. (9), the saturation concentration and 
the corresponding saturation pressure p s are related 
by: 


t 


E — 
evap " v 


( 4 ) 


The solution to Eq. (3) satisfying the initial condi- 
tion n( 0 ) - 0 is: 



(9a) 
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Let a circular void of radius R be suddenly created 
in the film and the dissolved gas start to diffuse into 
the void. It is intended to find how the pressure in 
the void increases with time. 

The variation of gas concentration in the solution 
is governed by the diffusion equation: 


i a_/ ac\ i 3C 
r 9r V dvj m D 9t ' 


R < r < ® 


( 10 ) 


The initial condition is: 


result, which will show that p is indeed a slowly 
varying function during the diffusion process. 

The mass flux into the void is given by Pick's 

law: 

’diff ■ D ®) r=R (18 

Using Eq. (17) to evaluate the concentration gradient, 
one finds: 


C(r,0) - C s 


(U) 


The boundary condition at large distance from the void 
is : 

C(® , t ) ♦ C s (12 ) 


With p ( t ) denoting the pressure in the void, Eq. (9) 
is applied as the boundary condition at r - R. let 
tdyn be a representative duration for the occurrence 
of cavitation, and ft the characteristic size of the 
region where diffusion effectively takes place. From 
Eq. (10) one deduces that 


a = ^Dtj- V /2 


dyn ) 


(13) 


For a gas-in-1 iquid system, typically D ~ 0(10’® 
m 2 /s); and in the operation of dynamically loaded bear- 
ings, typically t^yn ~ 0(10“^ s). Then Eq. (13) shows 
that ft * 0(10’® m). Hence, diffusion is effective 
only in a thin layer whose thickness is much smaller 
than R. Substituting 


r - R + ft?, 5 ~ 0(1) 

1 - *dyn T ’ T ~ 0(1) 


(14a, b) 


into Eq. (10) and retaining only the leading order 
terms, one obtains: 


3 2 C 

as 2 


ac 

3 t 


(15) 


The initial and boundary conditions are transformed 
into: 



The change* of the gas density in the void may then be 
expressed as: 


*R 2 h 2 ff Rhj diff (20) 

It is assumed that the gas in the void obeys the ideal 
gas law: 


p > (2i) 

Tn addition, a parameter t^.^ is defined to simplify 
the notation: 




V 1 


diff - D 


( 22 ) 


Substituting Eqs . (9a), (19), (21), and (22) into (20), 
one obtains the equation governing p: 


dp 

dt 


p s - p 

z ( l dif f 0 


1/2 


(23) 


The solution to this equation satisfying the initial 
condi tion p(0) ~ 0 is : 


EL 

p s 


■( t/t di f f) 2 


(24) 


At x ■ 0 , 
As £ ♦ ®, 

At l - 0, 



The parameter t ^ j f f signifies the time needed 
for the pressure in the void to reach 63 percent of 
p s , and may be taken as the characteristic time of dif- 
fusion. Its evaluation requires the knowledge of the 
coefficient of binary diffusion and Henry's constant. 
(16a,b,c) For aiT tn squalene, these quantities are found (Appen- 
dix) to be D =* 0.495x10“® m 2 /s and H - 13.4 MPa. 
Again by taking R =■ 10 mm and T - 293 K, one obtains 
from Eq. (22) that t rfiff - 0.324xl0 6 s. 

DISCUSSION 


The solution to this system can be readily found. In 
terms of the original variables, it is: 



This is an approximate solution because during its 
derivation p(t) is treated as constant. This quasi- 
steadiness approximation can be justified by the final 


The characteristic time of diffusion calculated in 
the above example appears unexpectedly long, and war- 
rents a comparison with the results of similar studies 
in the literature. Epstein and Plesset (1950) calcula- 
ted the time of completely dissolving a spherical air 
bubble in water. In their analysis the concentration 
of air at the bubble wall was assumed always saturated 
(C s ), and that in water (C| ) was treated as a parame- 
ter. The dissolving time was expressed as a function 
of the concentration ratio (Cj/C s ) and the initial 
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bubble radius (R 0 ) • In the case of Cj - 0, where 
diffusion was the fastest, the dissolving time was 
found to be 1.25 s for R 0 - 0.01 nan; 1.25xl0 2 s for 
R 0 - 0.1 tan; and 1.25xl0 4 s for R 0 - 1 nun. Recogni- 
zing from Eq. (22) that tjjff is proportional to the 
square of the bubble radius, one would expect a dis- 
solving time of 0(10® s) for a bubble radius of 10 mm. 
Krieger, Mulholland, and Dickey (1967) actually mea- 
sured the dissolving time. For oxygen bubbles of ini- 
tial radii in the neighborhood of 0.3 mm, they observed 
the dissolving time to vary from 6xl0 2 s to 0.45xi0 4 s, 
for C[ ranging between zero and 90 percent of satura- 
tion. Thus, the credibility of Eq. (22) is supported 
by these earlier results; they all point to the fact 
that diffusion is a very slow process. 

The load frequencies encountered in current oil 
film bearings fall in a certain range. The operating 
speed of the squeeze film dampers in gas turbine 
engines may be as high as a few ten-thousand revolu- 
tions per minute; while that of the main and connect- 
ing rod bearing in reciprocating internal combustion 
engines may be as low as a few hundred revolutions per 
minute. This range of operating speeds corresponds to 
a tjyn value from 1 ms to 0.1 s. Thus, t^yn is usu- 
ally much longer than the characteristic time of evapo- 
ration, and one may expect the presence of oil vapor in 
the cavitation bubble. On the other hand, t^yn is 
much too short compared with the characteristic time of 
diffusion. Hence, there is not enough time for the 
dissolved air to come out of the solution and find its 
way into the cavitation bubble. This analysis is con- 
sistent with the experimental results of Sun and Brewe 
(1990), where no difference in the cavitation pattern 
was observed between the degassed and the aerated oil 
throughout the cases studied. Since the diffusion of 
dissolved air does not play a role in the dynamic cavi- 
tation phenomenon, then, whether the oil has dissolved 
air in it or not should make no difference. 

It has often been reported* that, under dynamic 
bearing loads, persistent and numerous gas bubbles were 
found in the oil film, which eventually became a well- 
mixed two-phase fluid. The presence of these gas bub- 
bles has been attributed to both the liberation of 
dissolved gases in oil and the entrainment of ambient 
air. It seems clear that the former suspected cause 
can be ruled out by the result of the present analysis. 

CLOSURE 

Two formulas (Eqs. (4) and (22)) are derived which 
can be used conveniently to determine, respectively, 
the characteristic time of filling a void with the 
vapor of the surrounding liquid, and that of filling a 
void by diffusion of the dissolved gas in the liquid. 
These time scales can be compared with the characteris- 
tic time of bearing operation to determine the rele- 
vance of these processes with reference to the dynamic 
cavitation in oil film bearings. 

The computation shows that the evaporation process 
is usually fast enough to fill the cavitation bubble 
with oil vapor; whereas the diffusion process is much 
too slow for the dissolved air to liberate itself and 
enter the cavitation bubble. 
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APPENDIX 

CALCULATING D AND H FOR THE SOLUTION OF AIR IN 
SQUALENE 

Some needed properties of air and squalene are 
listed in Table 1. The coefficient of binary diffusion 
for gas in liquid may be calculated by using the Arnold 
equation (Simons and Ponter, 1975) which reads: 



where the empirical constant 0.010 is accurate at 
20 °C; the units of the quantities on the RHS of the 
equation are as shown in Table 1; and D is in cm 2 /s. 
Substituting the appropriate values of Table 1 into 
the equation, one obtains D = 0.495xl0“ 5 cm 2 /s. 

Henry’s constant can be estimated by the method of 
Osburn and Markovic (1969). The procedure is to first 
calculate a parameter “’parachor” for the liquid based 
on its molecular structure. From the parachor and the 
liquid’s molal volume, its surface tension can be eval- 
uated. Then from empirical correlations an Ostwald 
coefficient is obtained, which can be converted to 
Henry’s constant in unit of atmospheres. Following 
this procedure and with the use of Table 1, one obtains 
H - 132 atmospheres. 
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TABLE 1. - SOME PROPERTIES OF AIR AND SQUALENE 



Air 

Squalene 

Molecular weight 

Mg - 28.97 g/mole 

M t - 422.8 g/mole 

Mass of a molecule 

m - 70.2xl0~- 3 g 

Densi ty 


Pt » 0.81 g/ml 

Molal volume (at boil- 

Vg js 30 ml /mole 

V( = 522 ml /mole 

ing temperature) 
Absolute viscosity 

- 9-72 cP (at 20 °C) 

Parachor 


[P] - 1160.3 
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